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ABSTRACT 

We report the discovery of the 20.7 ms binary pulsar J 1952+2630, made using the distributed computing 
project Einstein® Home in Pulsar ALFA survey observations with the Arecibo telescope. Follow-up observa- 
tions with the Arecibo telescope confirm the binary nature of the system. We obtain a circular orbital solution 
with an orbital period of 9.4 hr, a projected orbital radius of 2.81t-s, and a mass function of / = 0.15 Mq by 
analysis of spin period measurements. No evidence of orbital eccentricity is apparent; we set a 2(7 upper limit 
e < 1.7 X lO^''. The orbital parameters suggest a massive white dwarf companion with a minimum mass of 
0.95 M0, assuming a pulsar mass of 1.4 Mq. Most likely, this pulsar belongs to the rare class of intermediate- 
mass binary pulsars. Future timing observations will aim to determine the parameters of this system further, 
measure relativistic effects, and elucidate the nature of the companion star. 

Subject headings: pulsars: general, pulsars: individual (J1952+2630), stars: neutron, white dwarfs 



1. INTRODUCTION 

' Albert-Einstein-Institut, Max-Planck-Institut fiir Gravitationsphysik, D- 
30167 Hannover, Germany 

^ Institut fur Gravitationsphysik, Leibniz Universitat Hannover, D-30167 
Hannover, Germany 

' Department of Physics, McGill University, Montreal, QC H3A2T8, 
Canada 

* Physics Department, University of Wisconsin - Milwaukee, Milwaukee, 
WI 53211, USA 

' Space Sciences Laboratory, University of California at Berkeley, Berke- 
ley, CA 94720, USA 

* Swinburne University, Center for Astrophysics and Supercomputing, 
Hawthorn, Victoria 3122, Australia 

' Astronomy Department, Cornell University, Ithaca, NY 14853, USA 

NAIC, Cornell University, Ithaca, NY 14853, USA 
' Columbia Astrophysics Laboratory, Columbia University, New York, 
NY 10027, USA 

'" Department of Physics and Astronomy, Franklin and Marshall College, 
Lancaster, PA 17604-3003, USA 

" Arecibo Observatory HC3 Box 53995, Arecibo, PR 00612, USA 

'^ Department of Astronomy and Radio Astronomy Laboratory, Univer- 
sity of California, Berkeley CA 94720, USA 

" Max-Planck-lnstitut fiir Radioastronomie, D-53121 Bonn, Germany 

'* Netherlands Institute for Radio Astronomy (ASTRON), Postbus 2, 7990 
AA Dwingeloo, The Netherlands 

" Astronomical Institute "Anton Pannekoek", University of Amsterdam, 
1098 SJ Amsterdam, The Netherlands 

'* Center for Gravitational Wave Astronomy, University Texas - 
Brownsville, TX 78520, USA 

Jodrell Bank Centre for Astrophysics, School of Physics and Astron- 
omy, University of Manchester, Manch., M13 9PL, UK 

Department of Physics, West Virginia University, Morgantown, WV 
26506, USA 

Also adjimct at the NRAO (National Radio Astronomy Observatory), 
Green Bank, WV 24944, USA 

Cardiif School of Physics and Astronomy, Cardiff University, Queens 
Buildings, The Parade, Cardiff, CF24 3AA, UK 

2' Department of Physics, Lafayette College, Easton, PA 18042, USA 

Albert-Einstein-Institut, Max-Planck-Institut fur Gravitationsphysik, D- 
14476 Golm, Germany 

NRAO (National Radio Astronomy Observatory), Charlottesville, VA 
22903, USA 

^ Department of Physics and Astronomy, University of British Columbia, 
Vancouver, EC V6T IZl, Canada 



Pulsars in short-period orbits with neutron stars or white 
dwarfs are invaluable tools for diverse areas of science. Pul- 
sars are precise clocks and enable very stringent tests of Ein- 
stein's theory of general relativity (e.g., Taylor & Weisberg 
1989 and Kramer & Wex 2009). These binary systems also 
provide unique opportunities to studying their properties as 
stellar-evolution endpoints (see review by Stairs 2004). The 
detection of relativistic effects like the Shapiro delay, most 
easily measured for highly inclined binary systems with a 
massive companion, can reveal the orbital geometry as well 
as the masses of the pulsar and its companion (Ferdman et al. 
2010). Precise mass measurements of the pulsar further our 
understanding of matter at (super)nuclear densities by provid- 
ing constraints on the possible equations of state, e.g., Demor- 
estetal. (2010). 

Here, we report the discovery and the orbital parameters of 
PSR J1952-I-2630. This pulsar is in an almost circular 9.4 hr 
orbit with a massive companion of at least 0.95 Mq, assuming 
a pulsar mass of 1.4 Mq. Thus, this new binary system is a 
good candidate for the measurement of relativistic effects like 
the Shapiro delay and could have an impact on all science 
areas listed above. 

2. DISCOVERY OF PSR J1952+2630 

PSR JI952H-2630 was discovered in Pulsar ALFA (PALFA) 
survey observations taken in 2005 August with the 305 m 
Arecibo telescope. The survey observations cover two sky 
regions close to the Galactic plane (\b\ < 5°). One region 
is in the inner Galaxy (32° ;$ ^ < 77°) and is observed 
with 268 s dwell time per observation; the other in the outer 
Galaxy (168° < £ < 214°) with 134 s long observations. 
The observations use ALFA, a cooled seven feed-horn, dual- 
polarization receiver at 1.4 GHz (Cordes et al. 2006). Signals 
are amplified, filtered, and down-converted. Then, autocorre- 
lation spectrometers, the Wideband Arecibo Pulsar Processors 
(WAPPs; Dowd et al. 2000), sum polarizations and generate 
spectra over 100 MHz of bandwidth with 256 channels every 
64 /is. 
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Data are archived at the Center for Advanced Computing 
(CAC), Cornell University. They are searched for isolated and 
binary pulsars in three independent pipelines: (1) a pipeline 
at the CAC searching for isolated pulsars and single pulses; 
(2) a pipeline using the PRESTO^^ software package oper- 
ating at several PALFA Consortium member sites, searching 
for isolated pulsars and binary pulsars with orbits longer than 
« 1 hr, and single pulses; and (3) Einstein® Home, search- 
ing for isolated or binary pulsars with orbits longer than 11 
minutes. 

Einstein @Home^^ is a distributed computing project. Its 
main goal is the detection of gravitational waves from un- 
known rapidly spinning neutron stars in data from the LIGO 
and VIRGO detectors (Abbott et al. 2009). Since 2009 March, 
about 35% of Einstein @Home compute cycles have been 
used to search for pulsars in radio data from the PALFA sur- 
vey. For Einstein® Home, volunteer members of the public 
sign up their home or office computers, which automatically 
download work units from project servers over the internet, 
carry out analyses when idle, and return results. These are 
automatically validated by comparison with results for the 
same work unit, produced by a different volunteer's com- 
puter. As of today, more than 280,000 individuals have con- 
tributed; each week about 100,000 different computers down- 
load work. The aggregate sustained computational power of 
0.38PFlops~^ is comparable to that of the world's largest 
supercomputers^^ . 

For Einstein@Home, raw data are transferred from the 
CAC to the Albert Einstein Institute, Hannover, Germany, 
via high-speed internet connections. There, servers dedis- 
perse the raw data into time series with 628 trial dispersion 
measure (DM) values up to 1002.4 pccm~^. For bandwidth 
and throughput reasons, the time resolution of the raw data 
is reduced by a factor of two to 128 /is. The dedispersed and 
downsampled time series are downloaded by the volunteers' 
computers over the internet and coherently searched for sig- 
nals from pulsars in circular orbits longer than 11 minutes. A 
detection statistic, the significance S = — log^g (p), is evalu- 
ated on a grid of parameter space points, where p is the false- 
alarm probability of the signal in Gaussian noise. A list of 
the 100 most significant candidates for each dedispersed time 
series is returned. After completion of all work units for a 
given observation, the results are post-processed on servers in 
Hannover, visually inspected and optimized using tools from 
the PRESTO software package and finally uploaded to a cen- 
tral database at Cornell. A more detailed account of the Ein- 
stein®Home pipeline is available in Allen et al. (2011 (in 
preparation). 

The 20.7-ms pulsar J1952H-2630 was found with a maxi- 
mum significance of 5 = 39.6 by visual inspection of the Ein- 
stein®Home results from a 268 s survey pointing with beam 
center at equatorial coordinates (J2000.0) a = 19''52™34.'5, 

5 — +26°31'14". The pulsar is detected most significantly 
at a DM of 315.4 pecm-^. The NE2001 model of Cordes 

6 Lazio (2002) with the given sky position implies a dis- 
tance of 9.41^ 4 kpc. The discovery observation exhibits a 
marginally significant, but large, barycentric period derivative 
P = 1-1(7) X IG^^ss"^ over the short observation time and 
indicates that PSR J 1952-1-2630 is in a short-period binary sys- 

http : / /www . cv . nrao . edu/~sransoin/presto/ 

http : //einstein .phys . uwm.edu/ 

http : //www. topSOO .org/list/2010/11/100 



tem. 

3. FOLLOW-UP OBSERVATIONS 

To obtain the orbital parameters of the binary system, 
follow-up observations with the Arecibo telescope were car- 
ried out between 2010 luly 29 and 2010 November 24. They 
were conducted mostly in coincidence with PALFA survey 
observations and used the central beam of the ALFA receiver 
with the Mock spectrometers^*^. These observations provide 
a total spectral range of 300 MHz in two overlapping bands 
of 172 MHz with 512 channels each, centered on 1.4 GHz at 
a time resolution of 65.476 /is. For our analysis, we use data 
from the upper band ranging from 1.364 GHz to 1.536 GHz, 
since the lower band tends to show more radio frequency in- 
terference from terrestrial sources. Most follow-up observa- 
tions cover an observation time Tobs ~ 600 s, though a few are 
of shorter duration. On 2010 July 29 and 2010 July 30 two 
longer (Tobs ~ 4200 s) follow-up observations with a time 
resolution of 142.857 /is covering the same frequency range 
in 2048 channels were carried out. The epochs of the follow- 
up observations range from MJDs 55407 to 55525. 

On 2010 August 19, gridding observations of the pulsar po- 
sition were performed at 2.1 GHz to improve the uncertainty 
of the discovery sky position. We used a square grid cen- 
tered on the discovery position with nine observation point- 
ings. The pulsar was found with equal signal-to-noise ratio 
in two of the gridding pointings. The improved sky position 
half-way between these two pointings is a = 19''52'"34.*4, 
5 = +26°30'14" with an imcertainty of w 1' given by the 
telescope beam size. 

Figure 1 shows the folded pulse profile of PSR J1952H-2630 
obtained from the 576 s follow-up observation on 2010 
September 19. The full width at half maximum duty cy- 
cle is 6%, corresponding to a width of the pulse of w^o = 
1.3 ms. The dispersive delay across a frequency channel is 
0.3 ms for the given DM, frequency resolution and central fre- 
quency. We calibrate the profile using the radiometer equa- 
tion to predict the observation system's noise level. Since 
the sky position is known to 1' accuracy, we can only de- 
rive limits on the flux density. We obtain a system equivalent 
flux density S'sys of 3.0-3.8 Jy (depending on the precise pul- 
sar position within the beam). The observation bandwidth is 
172 MHz, the observation time is 576 s, and the folded pro- 
file has N = 128 bins; then the expected off-pulse noise 
standard deviation used for calibration is 76-98 /iJy. The esti- 
mated period-averaged flux density of the pulsar at 1.4 GHz is 
70 /xJy < 51400 < 100 /iJy (depending on its position within 
the beam). 

4. SPIN PERIOD-BASED ORBITAL SOLUTION 

We obtain a solution for the orbital parameters based on 
measurements of the barycentric spin period and spin period 

derivative. The pulsar's orbital motion in a binary system 
causes changes of the apparent barycentric spin frequency 
over time due to the Doppler effect. For a circular orbit, the 
barycentric spin period P as a function of the orbital phase 
9 = r^orb {t — Tasc) measured from the time of passage of the 
ascending node Tasc is given by 

P (6) = Pa il + xn,,k cos (6)). (1) 

Here, Qgib = 27r/Porb is the orbital angular velocity for an or- 
bital period Porb, Po is the intrinsic barycentric spin period of 

http : / /www. naic . edu/~astro/mock . shtml 
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Figure 1. Folded pulse profile of PSR J1952+2630 at 1.4 GHz from the 
576 s follow-up observation on 2010 September 20. The resulting period- 
averaged flux density is 70 /iJy < S'1400 < 100 /ijy (depending on the 
precise pulsar position within the beam). The duty cycle is 6% of the pulsar 
period Pg. The dispersive delay across a frequency channel is 0.3 ms. 

the pulsar, and x = rsin (i) /c is the projected orbital radius 
in light-seconds. The radius of the pulsar orbit is denoted by 
r, i is the orbital inclination, and c is the speed of light. 

4.1. First estimates of the orbital parameters 

For each follow-up observation of length Tobs ^ 600 s, 
we obtain a pair of values, the barycentric spin period P 
and the spin period derivative P, using the PRESTO soft- 
ware package. For the two longer follow-up observations of 

Tobs ~ 4200 s, we obtain (P, P) pairs over seven contigu- 
ous stretches of approximately 600 s each. The acceleration 
is computed from a = cPP^^, yielding in total 34 (P, a) 
pairs. We compute the best-fitting orbital parameters by the 
method in Freire et al. (2001). We find x = 2.76(4) It-s and 
-forb = 9.3(1) hr; numbers in parentheses are estimated la 
errors in the last digit. 

4.2. Refining the orbital parameters 

The initial set of orbital parameters is extended and refined 
using the starting values from Section 4.1. We employ least- 
squares fitting of the full parameter set ilorb, Tasc, x and Pq 
with Equation (1). After convergence of the fit, the parameters 
are further refined and errors are estimated by Markov Chain 
Monte Cai-lo (MCMC) sampling (Gregory 2005) using an 
independence chain Metropolis-Hastings algorithm with flat 
proposal density functions in all parameters. Table 1 shows 
the refined orbital parameters for PSR J1952H-2630 obtained 
in this manner 

This simple model does not take into account the intrinsic 
spin-down of the pulsar Fitting for the spin down rate is not 
yet possible with the data reported here, because of the degen- 
eracy with an offset in sky position, which itself is not known 
to high enough accuracy. This effect results in a systematic 
error of the best-fit results at an unknown level. Future analy- 
sis of pulsar phase instead of spin period, will distinguish the 
spin-down and orbital effects. 

Figure 2 shows all measured spin periods P as a function 
of the orbital phase in the upper panel and the residuals in 
spin period in the lower panel. No clear structure is visible in 
the residuals and thus no striking evidence of non-zero eccen- 
tricity is apparent, justifying the choice of the circular orbit 



Table 1 

PSR J1952+2630 Parameters from a Spin-period based Analysis 

Parameter Value 

Right ascension, a (J2000.0) 19'"52."6 " 

DecUnation, 5 (J2000.0) +26° 30'" 

Galactic longitude, £ (deg) 63.25 " 

Galactic latitude, b (deg) -0.37 " 

Distance, d (kpc) 9A'^W 

Distance from the Galactic plane, 1 2 1 (kpc) 0.06 

Dispersion measure, DM (pc cm^'^) 315.4 

Period averaged flux density, S'1400 (/^Jy) < 100, > 70 

FWHM duty cycle (pulse width, w^q (ms)) 6% (1.3 ) 

Intrinsic barycentric spin period, Pq (ms) 20.732368(6) 

Projected orbital radius, x (It-s) 2.801(3) 

Orbital period, Po,b (d) 0.3918789(5) 

Time of ascending node passage, Tasc (MJD) 55406.91066(7) 

Orbital eccentricity, e < 1-7 X 10^^ (2o-) 

Mass function, /(Mq) 0.15360(1) 

Minimum companion mass, rric (Mq ) 0.945 

Median companion mass, med (Mq ) 1.16 

Note. — Thirty-four measurements of P and P between MJDs 55407 and 
55525 are used. The numbers in parentheses show the estimated Ict errors in 
the last digits. 

"The sky position was obtained from a gridding observation and has an 
accuracy of 1 ' given by the telescope beam size. 
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Figure 2. Top: points (with errorbars, too small to see) are the measured 
spin periods P as a function of the orbital phase computed from the orbital 
model in Table 1. The curve is the expected model spin frequency and the 
horizontal red line is at Po . Bottom: the residual difference between expected 
and measured spin frequency as a function of the orbital phase. No clear trend 
indicating non-zero eccentricity is evident. 

model. 

4.3. Upper limit on orbital eccentricity 

To set an upper limit on the eccentricity of PSR 
J1952h-2630's orbit, we use a spin-period model based on the 
ELLl timing model (Lange et al. 2001) which is suited for 
orbits with small eccentricity e. Then, Equation (1) in first 
order in e is modified to 

P [9] = Po [1 + xl7„-b (cos (61) + £1 sin (261) + cos (261))] , 

(2) 

where the additional constants are given by ei — e sin (a;) 
and £2 = ecos(a;), respectively. The angle lo is the longi- 
tude of the periastron measured with respect to the ascending 
node of the orbit. We now fit for the complete set of parame- 
ters (Po, X, fiorb, Tasc, £i, £2)- FoUowing the same method as 
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described in Section 4.2, we find the best fitting orbital solu- 
tion including the additional eccentricity parameters. MCMC 
sampling of the parameters space is also conducted for this 
spin period model (Equation (2)). 

The parameters (/qi 2;, riorbi Tlsc) only change slightly 
within the confidence regions obtained from the circular fit. 
The posterior probability distribution function (pdf) of the ec- 
centricity is consistent with a circular (e = 0) orbital model; 
from the posterior pdf we obtain a 2a upper limit on the ec- 
centricity of e < 1.7 X 10~^ 

5. DISCUSSION 
The mass function is defined by 

Att'^c^ (rric sin (i))"^ 

where nip and are the pulsar mass and the companion 
mass, respectively. Inserting the orbital parameters yields 
/ = 0.15360(1) M0 forPSR J1952H-2630. Assuming a pul- 
sar mass rup = 1.4 M©, we obtain a minimum companion 
mass mc > 0.945 M© for i = 90°. With i = 60° we obtain 
the median companion mass = 1.16 Mq. For a com- 

panion mass — 1.25 (smallest measured neutron star 
mass; Kramer et al. 2006), we obtain an inclination angle of 
i = 55.1°. 

No evidence of eccentricity e > 1.7 x 10"'^ is found from 
this analysis. A common envelope phase with mass trans- 
fer in the system's past can explain the almost circular orbit 
and the short orbital period. Also, the companion's progenitor 
most likely was not massive enough to undergo a supernova 
explosion; the supernova would have likely kicked the com- 
panion which almost guarantees a much higher orbital eccen- 
tricity. Thus, a neutron star companion is conceivable but un- 
likely. The companion most likely is a massive white dwarf. 
The high white dwarf companion mass, the compactness of 
the orbit, and the moderate spin period, indicate that the pul- 
sar most likely survived a common envelope phase (Ferdman 
et al. 2010). 

Low eccentricity, spin period, and high companion mass 
most likely place this system in the rare class of intermedi- 
ate mass binary pulsars (IMBPs; Camilo et al. 2001). The 
distance of PSR J1952H-2630 to the Galactic plane, \z\ w 
0.06 kpc, is comparable to that of the other five IMBPs (Ferd- 
man et al. 2010); this low scale height of the known IMBP 
population might be due to observational selection effects 
from deep surveys along the Galactic plane (Camilo et al. 
2001). 

Phinney (1992) derived a relation between the orbital pe- 
riod and eccentricity of low mass binary pulsars (LMBPs with 
companion mass 0.15 M0 ^ Wc ^ 0.4 M0). Figure 4 in 
Camilo et al. (2001) displays this relation for LMBPs and 
IMBPs. The LMBPs follow the theoretically predicted re- 
lation very well, while the IMBPs do not follow the same re- 
lation; they have higher eccentricities than LMBPs with the 
same orbital period. The slope however, seems to be very 
similar to the one for the LMBPs. This might suggest that 
there exists a similar relation for this class of binary pul- 
sars. An exact measurement of the orbital eccentricity of 
PSR J1952H-2630 from a coherent timing solution will add 
another data point that could help to test this hypothesis at 
short orbital periods. 

Detection of the Shapiro delay in PSR J1952H-2630 might 
allow precision mass estimates and strong constraints on the 
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Figure 3. Measurable Shapiro delay amplitude (not absorbed by Keplerian 
orbital fitting) for PSR J1952+2630 as a function of the pulsar mass and the 
companion mass. The dashed lines show constant inclination angles of 45°, 
60°, 75°, and 90° at the bottom. 

orbital geometry of the binary. Figure 3 shows the mea- 
surable Shapiro delay amplitude for PSR J1952H-2630 as a 
function of the pulsar mass and the companion mass; this 
is the part not absorbed by Keplerian orbital fitting (peak- 
to-peak amplitude of Equation (28) in Freire & Wex 2010). 
Preliminary timing observations using the Mock spectrome- 
ters at the Arecibo telescope have TOA uncertainties of ap- 
proximately 20 /is. These are currently unconstraining, be- 
cause no dedicated, deep observations at superior conjunction 
are available. Observations over a larger bandwidth and use 
of coherent dedispersion techniques with new instrumentation 
could improve TOA uncertainties further to w 10 /is. Thus, a 
detection of Shapiro delay requires relatively high inclination 
angles. If no Shapiro delay is detected, this will enable more 
stringent lower limits on the companion mass. These alone 
promise to be interesting, given the already known high min- 
imum companion mass. 

6. CONCLUSIONS AND FUTURE WORK 

We have presented a spin-period based analysis of the or- 
bital parameters of the newly discovered binary pulsar PSR 
J1952H-2630. The pulsar has an orbital period of 9.4 hr in an 
almost circular orbit with a projected radius of 2.81t-s. The 
mass function of f — 0.154 implies a minimum com- 
panion mass of 0.945 M0. Most likely, the companion is a 
massive white dwarf, although a neutron star companion can- 
not be excluded. 

The above described observations and further follow-up ob- 
servations will be used to derive a coherent timing solution for 
PSR J1952H-2630. This will provide a more precisely mea- 
sured sky position, orbital parameters, and values for the or- 
bital eccentricity and the intrinsic spin-down of the pulsar. 
This should enable a detailed description of this binary sys- 
tem and constrain its possible formation. A precise position 
would also enable searches for counterparts in X-ray, infrared, 
and optical wavelengths, although the large distance makes 
detections challenging. 

Furthermore, detection of Shapiro delay could be possible 
with further timing observations for high orbital inclinations. 
Given the already high minimum companion mass derived in 
this Letter, even a non-detection of the Shapiro delay could 
provide interesting, more stringent limits on the companion 
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mass and its nature. 

This pulsar is the second pulsar discovered by the global 
distributed volunteer computing project Einstein® Home 
(Knispel et al. 2010). This further demonstrates the value of 
volunteer computing for discoveries in astronomy and other 
data-driven science. 
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